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ABSTRACT 

General Relativistic (GR) Magnetohydrodynamic (MHD) simulations of black hole 
accretion find significant magnetic stresses near and inside the innermost stable cir- 
cular orbit (ISCO), suggesting that such flows could radiate in a manner noticeably 
different from the prediction of the standard model, which assumes that there are 
no stresses in that region. We provide estimates of how phenomenologically interest- 
ing parameters like the "radiation edge" , the innermost ring of the disc from which 
substantial thermal radiation escapes to infinity, may be altered by stresses near the 
ISCO. These estimates are based on data from a large number of three-dimensional 
GRMHD simulations combined with GR ray-tracing. For slowly spinning black holes 
(a/M < 0.9), the radiation edge lies well inside where the standard model predicts, 
particularly when the system is viewed at high inclination. For more rapidly spin- 
ning black holes, the contrast is smaller. At fixed total luminosity, the characteristic 
temperature of the accretion flow increases between a factor of 1.2 — 2.4 over that 
predicted by the standard model, whilst at fixed mass accretion rate, there is a cor- 
responding enhancement of the accretion luminosity which may be anywhere from 
tens of percent to order unity. When all these considerations are combined, we find 
that, for fixed black hole mass, luminosity, and inclination angle, our uncertainty in 
the characteristic temperature of the radiation reaching distant observers due to un- 
certainty in dissipation profile (around a factor of 3) is greater than the uncertainty 
due to a complete lack of knowledge of the black hole's spin (around a factor of 2) 
and furthermore that spin estimates based on the stress-free inner boundary condition 
provide an upper limit to a/M . 
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1 INTRODUCTION 

Recent years liave seen a rapid growth in the specificity 
and detail with which we attempt to describe accretion onto 
black holes. In the early days of the subject three decades 
ago, we were c ontent with the simplest o f models (the 'stan- 
dard' model, iNovikov fc Thornel [l973l ') , hereafter NT, in 
which the accretion flow was assumed to be time-steady, ax- 
isymmetric, geometrically thin, and following circular Ke- 
plerian orbits at all radii outside the innermost stable cir- 
cular orbit (the ISCO, located at radius rms). Relying on 
heuri stic arguments framed in a purely hydrodynamic con- 
text (|Page fc Thornelll974l ). it was further assumed that all 
stresses ceased inside the ISCO, so that the inner edge of the 
disc could be described as falling precisely at that radius. 
Dimensional analysis was the basis for any link between the 
inter-ring s tresses essential to accret ion and local physical 
conditions (jShakura fc SunvaevlllQT^ '). 



Today, there is tremendous effort to obtain and 
interpret direct observational diagnostics of the inner- 
most parts of accretion flows onto black holes. Appar- 
ently relativistically-broadened Fe Kq proflles can be dis - 
cerned in numerous examples (|Revnolds fc NowakI |2003| ) . 
Extensive efforts are made to fit det ailed disc spectral 
models to observed continuum spectra |Gierlihski fc Done 



2004lDavis et al.ll2005l:IShafee et al.ll2006l : lMcClintock et al 



2009 : iHui fc Krolikl 120071 ). These spectral fits can be 
used to infer the mass of the central black hole; 
both methods may be used to const r ain the spin of 
the b lack hole (iMakishima et all |200G|: iMiller fc ColbertI 



20041: iMiniutti et all 12004 iBrenneman fc Revnoldsl l20o'i 



Shafee et all l2006l : iMcClintock et al.l bood ) . Some hope to 
use relativistic fluid dynamics to deflne normal modes of disc 
oscillation that could then also be used to constrain the cen- 
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tral black hole's spin (jWagoner et alj 1200 ll : iRezzolla et al.l 
l2003l : lKatd[2004l ). 

At the same time, there has also been much progress 
on the theoretical side. A strong consensus now supports 
the idea that accretion stresses are the result of correlated 
MHD turb ulence, driven by a perv asive magneto-rotational 
instability (|Balbus fc Hawlevlll998l '). Numerical simulations 
developing this idea have given us a detailed and quanti- 
tative description of the vertical profiles of pressure and 
density inside disc s (|Miller fc Stonell2000l : lHirose et al.|[200^ : 
iKrolik et aLll2007l ). as well as provided us with detailed pic- 
tures of how their p roperties vary smoothly across the ISCO 
l|Krolik et al. I [20051 ). Thes e simulations h ave vindicated the 
prescient remark made by lThorn3 (|l974 ): 
"In the words of my referee, James M. Bardeen (which echo 
verbal warnings that I have received from Ya. B. Zel'dovich 
and V.F. Schwartzman), 'It seems quite possible that mag- 
netic stresses could cause large deviations from circular or- 
bits in the very inner part of the accretion disk and change 
the energy-angular-momentum balance of the accreting mat- 
ter by an amount of order unity'." 

Indeed, when the black hole rotates, significant mag- 
netic stresses can be found throughout the ac cretion fiow, 
all th e way to the edge of the event horizon (jKrolik et al.l 
l2005l ). 

Work on specific dynamical pictures of accretion has 
stimulated a reexamination of the simple picture that discs 
have sharp edges, within which little of interest happens. 
Although it is true that there are qualitative changes across 
the ISCO, they do n o t hap pen discontinuously. As argued 
by iKrohk fc HawWl iooi), what one means by "edge" 
depends on the question asked. For example, the "refiec- 
tion edge", the edge outside of which most of the ob- 
served Fe Ka and Compton reflection photons are created, 



ISCO dRevnolds fc Begelman 


1 19971: Krohk fc Hawlev 20021: 


iBrenneman fc Reynolds! 200C 


; Reynolds fc Fabian 200^)— 



its exact position depends on the density and optical depth 
of the gas in that region, and on the intensity of ionizing ra- 
diation striking it. Similarly, the "radiation edge" (the edge 
inside of which little of the total luminosity emitted by the 
flow escapes to infinity) should be near the ISCO, but is not 
necessarily identical to it. Its position depends both on the 
profile of dissipation and on the ability of photons to escape 
to infinity, and to do so without excessive loss of energy to 
gravitational redshift. 

In this paper we examine the infiuence of magnetic 
stresses at and inside the ISCO on the apparent size, lu- 
minosity and characteristic temperature of black hole accre- 
tion discs. This effort is important to black hole phenomenol- 
ogy because so many observational diagnostics depend upon 
these three parameters. A prime example is provided by at- 
tempts to use spectral fits to constrain black hole mass and 
spin. This program rests upon the idea that the observed 
luminosity, L, is essentially thermal, so that it may be char- 
acterized by an effective radiating area and a characteristic 
temperature; L = A{a/ M,9)T*f^^^. Because both T^har and 
L can be measured, A can be inferred through this relation. 
With a theoretically-supported connection between A and 
a/M (and some other constraint on 9), the inference of A 
leads to an inference of a/M . C urrent efi'orts connec t Tchar 
to M, M, and a/M using the iNovikov fc Thornd (119731 ) 



model for the radial dissipation profile. This model depends 
in an essential way on the assumption that all stresses cease 
at and within the ISCO, whose location (at Boyer-Lindquist 
radial coordinate r/M = rms) is a function only of the black 
hole mass and spin. This temperature is further adjusted 
by an appropriate correction for gravitational and Doppler 
energy shifts and a "col or temperature corre ction" due to 
opacity effects (see e.g. iDone fc David l2008h . The appar- 
ent size of the disc. A, is likewise found from the dissipa- 
tion profile of the Nov ikov-Thorne model. As first noted by 
IPage fc Thornd (|1974| ). significant magnetic forces undercut 
the rationale for the zero-stress boundary condition; the sim- 
ulation data we discuss here shows quantitatively how these 
magnetic forces alter the connections between both A and 
Tchar and the black hole spin parameter. Although more 
work is needed to explore fully these new effects, in this pa- 
per we begin the discussion of how they can influence these 
inferences. 

The simulations reported here employ full general rel- 
ativity and three-dimensional MHD, so that whilst their 
treatment of angular momentum flow and inflow dynamics 
is quite accurate, they do not directly track dissipation. In 
an accretion disc, energy is extracted from orbital motion 
and transformed into kinetic and magnetic energy on the 
largest scales of the turbulence. Subsequently, this energy 
cascades down to a dissipation scale (either viscous or resis- 
tive) where it is finally thermalized. Current simulations can 
describe well the first stages of this process, but can mimic 
only indirectly the last step: grid- level effects intervene at 
lengthscales far larger than the physical scale of dissipation. 
In fact, the simulation code whose data we will use solves 
only the internal energy equation, and makes no attempt to 
follow dissipative energy losses except those associated with 
shocks. 

We proceed by instead making a plausible ansatz for 
heating within the disc that can be determined a posteriori 
from simulation data. As we shall see, our results depend 
primarily on the qualitative fact that dissipation continues 
smoothly across the ISCO, and on the nature of photon 
trajectories deep in a relativistic potential; for this reason, 
we believe that a non-rigorous, but physically-motivated, 
ansatz will not be misleading. We opt for the simplest choice: 
a connection between the heating rate and the stress that 
follows from the s tandard model for energy conservation in 
an accretion disc JPage fc Thorndll974l : iBalbus et al1ll994l : 
iHubenv fc Hubenvlll998D ~ 

We can also investigate the importance of enhanced 
stress in a semi-analytic fashion that bypasses most of 
the limitations of the current simula t ions. We employ the 
model formulated by lAgol fc KroTi^ (|2000l ) (hereafter the 
AK model) that allows for non-zero stress at the ISCO but 
otherwise computes the dissipation profile using the same 
approach as the standard model. The AK model cannot be 
extrapolated into the plunging region, and has in common 
with the NT model a fixed disc size. Its unique feature is the 
enhanced total dissipation due to the nonzero stress at the 
ISCO; simulation data provides the single parameter needed 
to calibrate the model. The AK model, therefore, provides 
an important link between the standard model and the full 
simulation results and allows us to gauge the appropriate- 
ness of the radiation ansatz employed for the latter. 

To relate dissipation rates to radiation received at infin- 
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Table 1. Simulation Parameters 



Name a/M r+/M nn/M r^s/M Field TavgXW^M Originally Presented in 



KDOb 


0.0 


2.00 


2 


104 


6.00 


Dipole 


6-8 


De Villiers et al. ('2003') 


KDOc 


0.0 


2.00 


2 


104 


6.00 


Dipole 


8-10 


Hawlev & Krolik (2006) 


QDOd 


0.0 


2.00 


2 


104 


6.00 


Quadrupole 


8-10 


This work 


KDIa 


0.5 


1.86 


1 


904 


4.23 


Dipole 


6-8 


De Villiers et al. ('2003') 


KDIb 


0.5 


1.86 


1 


904 


4.23 


Dipole 


8-10 


Hawlev & Krolik (2006) 


KDPd 


0.9 


1.44 


1 


503 


2.32 


Dipole 


6-10 


De Villiers et al. (2003) 


KDPg 


0.9 


1.44 


1 


503 


2.32 


Dipole 


8-10 


Hawlev & Krolik (2006) 


QDPa 


0.9 


1.44 


1 


503 


2.32 


Quadrupole 


8-10 


Beckwith et al. (2008) 


TDPa 


0.9 


1.44 


1 


503 


2.32 


Toroidal 


20-22 


Beckwith et al. (2008) 


KDG 


0.93 


1.37 


1 


458 


2.10 


Dipole 


8-10 


Hawlev & Krolik (2006) 


KDH 


0.95 


1.31 


1 


403 


1.94 


Dipole 


8-10 


Hawlev & Krolik (2006) 


KDJd 


0.99 


1.14 


1 


203 


1.45 


Dipole 


8-10 


This work 


KDEa 


0.998 


1.084 


1 


175 


1.235 


Dipole 


6-8 


De Villiers et al. (2003) 


KDEb 


0.998 


1.084 


1 


175 


1.235 


Dipole 


8-10 


This work 


QDEb 


0.998 


1.084 


1 


175 


1.235 


Quadrupole 


8-10 


This work 



Here a/M is the spin parameter of the black hole, is the horizon radius, rin is the innermost radius in the computational grid, field 
is the initial field topology in the torus and Tavg x 10'^ is the time-interval over which simulation data was averaged. For reference, we 

note where individual simulations were originally presented. 



ity, we must perform one additional calculation using three 
further assumptions: that the dissipated heat is efficiently 
converted to photons, that these photons emerge from the 
accretion flow very near where they are created, and that 
they are radiated isotropically in the fluid frame. The flrst 
two assumptions are equivalent to requiring the timescales 
for dissipation, radiation, and photon diffusion to be shorter 
than the inflow timescale for all fluid elements. The third, 
while not strictly justifled, is the simplest guess we can make. 
Given those assumptions, we use a general relativistic ray- 
tracing code to relate the luminosity radiated by each fluid 
element to the luminosity received by observers located at 
different polar angles far from the black hole. 

A further result of this calculation is a new estimate of 
the radiative efficiency of accretion. The traditional calcu- 
lation of this quantity follows directly from a primary as- 
sumption that there are no forces inside the ISCO and two 
additional assumptions, that the radiation is prompt and all 
of it reaches infinity. We improve upon this traditional esti- 
mate in two ways: we allow for dissipation associated with 
the stresses we measure at and inside the ISCO; and we 
calculate the radiated energy (even within the NT model) 
that actually reaches distant observers. However, we do not 
regard our result as sufficiently final or complete to give it 
much weight, as it is likely to be more model- and parameter- 
dependent than our placement of the radiation edge. One 
reason for downplaying this result is that we find it neces- 
sary to omit any estimate of dissipation outside the main 
part of the accretion flow. We have less confidence that our 
dissipation prescription is appropriate in the jet, or even in 
the disc corona, than we do when it is applied to the main 
disc body and plunging region. Moreover, radiation from 
these lower-density regions is less likely to be thermalised 
and contribute to the radiation usually identified with the 
disc continuum. 

The rest of this paper is structured as follows. In ^ 
we briefly review the numerical scheme employed to solve 
the equations of GRMHD in the simulations whose data we 
use, give an overview of the parameters of the simulations 
included in this work, and describe our general relativis- 



tic ray-tracing code. In ^ we contrast the dissipation rate 
distributions predicted by the standard model, its Agol & 
Krolik extension, and our simulation-based ansatz. In ^we 
compute the luminosity at inflnity predicted by each of these 
dissipation proflles and discuss their consequences for the lo- 
cation of the radiation edge and the characteristic temper- 
ature of the accretion flow. Finally, in ^ we draw specific 
conclusions from our results and describe their implications 
for black hole phenomenology. 



2 NUMERICAL DETAILS 
2.1 GRMHD Simulations 

The calculations presented here continue the analy- 
sis of a progr am of black ho l e accr etion disc simula- 
tions begun in De Villiers et al. (|2003|) and continued by 



iHirose et al] (l2004l. iDe Villiers et al.1 (120051), iKrolik et al] 



(|2005l ). lHawlev fc Krolilj (|2006D and lBeckwith et al.l (120081) 
The simulation code is described in De Villiers fc HawlevI 



(|2003al ) . This code solves the equations of ideal non-radiative 
MHD in the static Kerr metric of a rotating black hole us- 
ing Boyer-Lindquist coordinates. Values are expressed in 
gravitational units {G = M — c = 1) with line element 
ds^ = gttdt'^ + 2gt^dtd(j) -\- Qrrdr^ + geedd^ + g^dydij}'^ and 
signature (— , -|-, -|-). The determinant of the 4-metric is 
Q^V7) where a = (— s")"^^'^ is the lapse function and 7 is 
the determinant of the spatial 3-metric. 

The relativistic fluid at each grid zone is described 
by its density p, speciflc internal energy e, 4- velocity u^, 
and isotropic pressure P. The relativistic enthalpy is h = 
1 -I- e -I- P/p. The pressure is related to p and e via the equa- 
tion of state for an ideal gas, P = pe{T — 1). The mag- 
netic fleld is described by two sets of variables. The first 
is the constrained transport magnetic field B' = [ijk]Fjk, 
where [ijk] is the completely anti-symmetric symbol, and 
Fjk. are the spatial components of the electromagnetic field 
strength tensor. From these are derived the magnetic field 
four-vector, (47r)^/^6'' = "F^^u,,, and the magnetic field 
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scalar, = h^h^. The electromagnetic component of the 
stress-energy tensor is T^^^^ = ig*^" 1 1 6| | ^ + uf'u" 1 1 b| | ^ - . 

The initial conditions for the simulations consist of an 
isolated hydrostatic gas torus orbiting near the black hole, 
with an inner edge at r = 15Af, a pressure maximum located 
at r « 25M, and a (slightly) sub-Keplerian distribution of 
angular momentum throughout. Parameters for all of the 
simulations analyzed in this work are summarised in Table 
[H We use three different initial field configurations. Models 
labeled KD have an initial dipole field that lies along sur- 
faces of constant pressure within the torus. The initial field 
in the QD models is a quadrupolar configuration consist- 
ing of field of opposite polarity located above and below the 
equatorial plane in the torus. The TD model's initial field 
is purely toroidal. Many of these simulations are studied in 
more detail in the referenced work given in the last column 
of the table. A few simulations are new to this work and 
these were computed with the latest vers ion of the GRMHD 
code as described bv iDe VilliersI (|2006t l. Whenever we re- 
quire time-averaged data, we average over the times Tavg 
shown in Table[Tl from which we have full three-dimensional 
snapshots taken every 80M. Thus, 26 time-samples go into 
each time-average. 

Each simulation uses 192 x 192 x 64 (r, 6, (j)) grid zones. 
The radial grid extends from an inner boundary located just 
outside the black hole event horizon (see Table [l]for the pre- 
cise location in each case), to the outer boundary located at 
rout = 120Af in all cases. The radial grid is graded accord- 
ing to a hyperbolic cosine function in order to concentrate 
grid zones close to the inner boundary, except in the TDPa 
model where a logarithmic grid was used. An outflow condi- 
tion is applied at both the inner and outer radial boundary. 
The 9-grid spans the range 0.0457r ^ 0.9557r using an 
exponential distribution that concentrates zones near the 
equator. A reflecting boundary condition is used along the 
conical cutout surrounding the coordinate axis. Further dis- 
cussion of the r, 6 gridin g and boundary conditions can be 
found in lHawlev fc Kroli k (.2006,). Finally, the (ji-grid spans 
the quarter plane, 0^0^ 7r/2, with periodic boundary 
conditions in <j). The use of this restricted angular domain 
significantly reduces the computational requirements of the 
simulation ( for further discussion of the effects of this re- 
striction see iDe Villiers fc Hawlevll2003bl) . Most of the sim- 
ulations were run to either 8000 or lO^M, which corresponds 
to approximately 10 or 12 orbits at the initial pressure max- 
imum. The toroidal simulation was run to time 29000M due 
to the longer timescale take n for accretion to commence in 
this case (see iHawlev fc Krolik 2002,) . For each simulation 
the time step At is determined by the extremal light cross- 
ing time for a zone on the spatial grid and remains con stant 
for the entire simulation (|De Villiers fc Hawlevl[2003al ). 

In all the models the initial evolution is driven by the 
MRI and by growing magnetic pressure due to shear am- 
plification of poloidal field (if present) into toroidal. By the 
end of the simulation, a quasi-steady state accretion disc ex- 
tends from the hole out to r ~ 20M. Beyond this radius the 
net mass motion shifts to outward flow as it absorbs angu- 
lar momentum from the inner disc. In this paper, therefore, 
we focus our attention on the region inside 20M and on 
late times after the quasi-steady disc has been established. 
Whenever we require time-averaged data, we average over 
the times Tavg shown in Table [TJ from which we have full 



3-d snapshots taken every SOAf. Thus, 26 time-samples go 
into each time-average. 

2.2 Ray-Tracing 

Transformation into the reference frame of a distant observer 
was accomplished by means of a ray-tracing code. We assign 
to this observer a coordinate system ( "photographic plate" ) 
defined by the photon impact parameters (a, f3), which can 
be s imply related to the photon's constants of motion, (A, q) 
via (|Bardeen et al.lll972h : 

a = P = ±^q + a'^cos^Qo ~ X^cot^Qo (1) 

smfc^o 

where Qo is the co-latitude (inclination) of the distant 
observer. Once (A, q) are known, then the photon 4- 
momentum, at the observer can be constructed. Photon 
paths were then traced from a distant observer to the sur- 
face of the accretion disc by integration of the null-geodesic 
equation; 



\- Y \ — 

da^ '"^ ds da 



(2) 



where is the spacetime coordinate of the photon, s is 
an affine parameter, FJ^^ are the connection coefficients and 
the dx^ /da are specified initially by the at the distant 
observer. Equation [2] may be recast as a set of coupled 
first-order differential equations. We solve this set by ap- 
plying the fifth-order Runge -Kutta integrator described by 
iBrankin fc Shampinj (|l991h as implemented in the NAG 
FORTRAN Library (Mark 21). The local error in the in- 
tegration was kept to 10^*. Where the geodesic passed 
through a radial turning point, the location of the turning 
point found by the integrator was com pared with its ana- 
lytic value (see iBeckwith fc DondlioOSl ) and was found to 
agree to the order of the local error in the integration. The 
integration of the geodesic was terminated when the pho- 
ton intersected the disc surface, which was located in the 
9 = n /2 plane. Once the intersection with the disc surface, 
a;^^j.j(a, /3) is known, the flux measured by the distant ob- 
server can be calculated. By treating the emission as a line 
flux, a given radius within the disc can be associated with 
an observed flux by 



F{a,l3) = /[<„,f(a,/3)]Q[<^,f(a,/3)]5(S, - gE, 



(3) 



The term g is the Doppler and gravitational energy shift 
of the photon between the disc and the distant observer, 
found by projecting the photon four-momentum onto the 
tetrad describing the fluid rest frame (see Appendix A). For 
the standard disc models the fluid motion is assumed to 
be purely Keplerian, that is = l/{a + r^^^), whilst for 
fluxes calculated from the simulations the fluid motion on 
the emission surface is obtained directly from the data. 
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3 DISSIPATION IN THE FLUID FRAME 

The radial dependence of the fluid-frame dissipation rate per 
unit disc surface area in the standard model is 



3GMM 

4yj-j.3 



RB.{r), 



(4) 



where RR{r) is a function encapsulating all of the relativistic 
effects relevant for disc dynamics and the effect of the as- 
sumed stress-free inner boundary condition at the ISCO. In 
the Newtonian limit for a Keplerian potential, RR{r) takes 
the familiar form: 



RR{r) = 1 - (r^/r) 



1/2 



(5) 



10.00 



1.00 



0.10 



0.01 




0.001 



0.010 



0.100 



1.000 



1-a/M 



where rm is the radius in the disc at which the stress is zero. 
Q is defined in the frame that comoves with the surfaces of 
the disc (the "disc frame"). Determining the form of this 
function as measured by a distant observer then requires a 
transformation into that observer's reference frame, accom- 
plished by explici tly tracing photons from the emitter to the 
observer (see e.g. iBeckwith fc Done|[2004l ). 

The NT dissipation profile is derived from conserva- 
tion laws applied in a time-steady context. We require 
the instantaneous dissipation in a non-stationary flow. To 
estimate it, we begin with the vertically-integrated and 
azimuthally-averaged fluid-frame Maxwell stress VV^'^j. Fol- 
lowing iKrohkeFal] (120051 ) ■ we find the fluid-frame Maxwell 
stress from the following expression; 



/dis 



rpfj. 



-M =(0)-'^ (EM) 



/ dx('l''>{r,e = -k/2) 



(6) 



where e'f , are the basis vectors describing the rest frame of 

(") 

the fluid (see Appendix A) and da;'^' = e\!^^ dx" is the fluid- 
frame co-ordinate element. The subscript "disc" denotes 
that we are including only contributions to the integral from 
bound material {—hut < 1) that lies within one density- 
scale height of the midplane (we deflne the density scale 
height 9h at a given r and (j) by p{r,(f>,9h) = p(r, 0, 7r/2)/e). 
This form of the stress facilitates comparison with the con- 
ventional representation of the ver tically-integrated stress 
given by iNovikov fc Thornd (|l973h : we discuss later why 
we restrict the integration to matter within a single density 
scale height of the plane. Plots showing the radial proflle 
of Wf^j com pared to the conventional stress are given in 
iKrolik et all (|2005| ) for four different black hole spins and 



Beckwith et al.l (|2008l ) for several different magnetic fleld 



topologies. In every case examined in both these studies, the 
Maxwell stress in the fluid frame extends all the way to the 
event horizon whenever the black hole rotates and reaches 
deep into the plunging region even when a/M = 0. 

To determine the observational implications of this ad- 
ditional stress we must compute its associated dissipation. 
As discussed above, this requires adopting a model for dis- 
sipation; we consider two. The first is the AK model, which 
extends the standard disc approach to allow for a non-zero 
stress at the ISCO. The second is a model, derived below, 
that determines the local dissipation in terms of the local 
stress and gradients of the fiuid velocity. 



Figure 1. Time-averaged fractional increase in total dissipation, 
(ADak)i derived using the Agol-Krolik model and the simulation 
Maxwell stress at the ISCO. Error bars show ±1 standard devia- 
tion about the average and reflect the intrinsic variability in the 
stress at the ISCO se en in the simulation. Th e dipolar field sim- 
ulations described by Id^ Villiers et al are shown by star 
symbols, those described by iH^leTfc Kro lik (2006) are marked 
by diamonds and the new high spin cases computed for this paper 
are marked by triangles. Simulations with quadrupolar topologies 
are denoted by squares. The toroidal field simulation (a/M = 0.9) 
is denoted by a cross. Simulation data was time-averaged over the 
intervals given in Table [T] 



3.1 Dissipation Derived from the Agol &l Krolik 
Model 

AK showed how the inner boundary condition for the stan- 
dard NT model can be modified to allow for arbitrary stress 
at the ISCO. Parameterizing the stress at the ISCO in terms 
of the resultant increase in radiative efficiency Ae yields the 
following expression for the vertically- integrated r — <j) stress 
in the fluid frame: 



dzS, 



(r)(0) 



MP. 
1^ 



,3/2 



+ RT{r) 



(7) 



where S(j.){<f>) is the r — (j> component of the fluid-frame vis- 
cous stress tensor; C(r), D{r) are relativistic correction fac- 
tors and Rrir) is similar in nature to RR{r), reducing to 
an identical expression in the Newtonian limit. The corre- 
sponding expression for the fluid-frame dissipation rate Qak 
is: 



3GMM 

j,3 



3/2 
' ms 



C{r) 



r.1/2 



+ RR(r) 



(8) 



This semi-analytic prescription lacks only a specification of 
the increase in radiative efficiency. This is supplied from the 
simulation data by setting — J dzS(^r){^){z) = ^(0) 
ISCO for each time-step in the data-set and then inverting 
Equation [7] to determine Ae. Lastly, we average Ae over 
the different time-steps. It is important to recognise that 
the nature of the AK model is such that stress data from 
one location alone — the ISCO — are enough to determine the 
stress (and dissipation) everywhere else because all locations 
in the disc outside Vms are linked through the assumptions 
of time-steadiness and axisymmetry. 
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It is also important to recognise that there are certain 
intrinsic limitations built into the AK model, inherited from 
its NT parent. These are best illustrated by returning to 
its basis, the relativistic conservation of momentum-energy, 
expressed by the requirement that VyT'^" = 0. In both of 
these models, the stress tensor was defined to have the form 



■ phu^'u-' + Pg^"" + 5""' + u^'g" + u^q" 



(9) 



where 5'"' is the part of the stress tensor responsible for 
inter-ring torques and is the energy flux four-vector. En- 
ergy flux contributes to the stress tensor in proportion to 
the four-velocity because energy can be directly advected 
with the material. It is assumed that both 5''"^ and q'^ are 
orthogonal to the four velocity: Uf^S'^" — and tt^g^ — 0. 
In the case of the stress tensor, this choice was made with 
classical viscosity in mind and also ensures that heat advec- 
tion can be cleanly separated from other kinds of stress. In 
the case of the energy flux four-vector, orthogonality with 
the four- velocity ensures true stationarity. Lastly, 5'"' is also 
assumed to be symmetric. In the context of MHD stresses, 
we have 



^l^^u-" + {\\bf /2)g'"' - bf'b'' , (10) 
which is manifestly symmetric, but is not orthogonal to u^^: 



^ (EM) 



'^M-' (EM) 



(lifeilV2K- 



(11) 



In other words. Maxwell stresses fail to be orthogonal to the 
fluid four-velocity to the degree that field energy is carried 
by fluid motions. 

Moreover, both the AK and NT models also assume 
that Ur = U0 = Q in the (Boyer-Lindquist) coordinate frame. 
Combined with the orthogonality between the four-velocity 
and the stress tensor, this assumption places rather spe- 
cial conditions on the relation between 6'* and fe'. Although 
Ur and ue are very small compared to and ut in the 
disc body, this assumption begins to fail as the ISCO is 
approached, and becomes entirely invalid in the plunging 
region. 

In any event, independent of the method used to find 
it, Q{r) can be used to compute the instantaneous total 
dissipation by integrating 



D = 



(12) 



where dx''^^ = ei dx'^ are the fluid frame differential coor- 
dinate elements in the 6 = Tr/2 plane. The time-averaged 
fractional increase in total energy dissipated due to the ad- 
ditional stress at the ISCO is then given by (ADi): 



(AA) 



(13) 



Here the index i could be either AK or MW, and all 26 
of our full three-dimensional snapshots are used for these 
time-averages. 

Making use of the data from all our simulations, we 
plot the fractional increases in dissipation for the AK model 
in Figure [T] This model predicts that nonzero stress at the 
ISCO produces a dissipation rate that is ~ 10% greater than 
in the NT model for zero spin, and increases steadily as the 
spin increases; when a/A4 — 0.998 the total dissipation rate 
is several times greater than in the NT model, with the 
specific factor depending on the initial magnetic topology. 



Generally speaking, non-dipolar magnetic field topologies 
lead to fractional increases that are smaller than their dipo- 
lar counterparts. This contrast arises from a distinguishing 
characteristic of the dipole simulations for high spin holes: 
in the region near r^s, the stress rises more rapidly toward 
small radius for this field ge ometry than for quadr upolar or 
toroidal fields (see Fig. 4 in lBeckwith et all (120081 ) ). 

3.2 Dissipation Derived from the Local 
Stress-Energy Tensor 

To derive dissipation rates directly from instantaneous and 
local simulation data, we return to the momentum-energy 
conservation equation V^T'^" = 0. If we define T*^" in the 
way it was defined for the NT and AK models, projecting 
the momentum-energy conservation equation onto ?t,i yields: 



+^i^q" + q'^u'V^u^ = 



(14) 



The GRMHD code evolves the internal energy equation 

V^(pe^t'')+PV^u^ =0, (15) 

without either a viscous stress (except for artificial bulk vis- 
cosity) or a radiative energy flux. We can then use the or- 
thogonality of q'^ and , along with the ansatz that energy 
liberated by the stress is promptly radiated, to write the 
energy balance equation as 



S^'V.ii^ -f V.g" = 0. 



(16) 



This description of the dissipation is consistent with the 
spirit of the standard disc approximation that the stress 
transporting angular momentum is the same str ess that lo- 
cally heats the disc (jBalbus fc Papaloizoul [l999[ l . However, 
if S*^^ is calculated from the electromagnetic stress-energy 
tensor, T^^y^y it is not self-consistent for the reasons dis- 
cussed in the context of Qak- 

To solve this equation, we assume that the energy flux 
escapes perpendicular to the equatorial plane so that g** = 
(0,0,0, g*) in the coordinate frame. We also assume that 
the flow is time-steady and axisymmetric (consistent with 
our use of time- and toroidally-averaged values) and that 
= in the coordinate frame. With these assumptions, we 
can expand the energy balance equation to yield: 



de{a^/jq ) = - [S ''dr{ct^ut) 

where we have used the identity Vp(/u^) 
(Q^/7)~"'^9,i(ay^/«^). Solving for yields: 



(17) 



, [5'*'■9r(a^/7Wt) 

Jdisc 



(18) 



where again the subscript "disc" denotes that we include 
only contributions to the integral from bound material 
{—hut < 1) that lies within one density-scale height of the 
midplane. The fluid frame dissipation rate computed from 
the electromagnetic stress-energy tensor , Qmw, is found 
from Qmw = e]^^ q^ , which for u'^ = (it', u*", 0, u**) yields: 



(19) 
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KDOc; a/M = 0.0 



KDIa; a/M = 0.5 




KDPg; a/M = 0.9 




KDEb: a/M = 0.998 



Figure 2. Time-averaged radial profile of Qmw (solid lines) compared to Qnt (dashed lines) and Qak (dot-dash lines) for (from left 
to right) dipole simulations KDOc {a/M = 0, KDIb (a/M = 0.5), KDPg {a/M = 0.9) and KDEb {a/M = 0.998). Dashed Unes denoted 
±1 std. deviation from the mean Qmw- Simulation data was timc-averagcd over the periods given in Table[T] 



KDPg; a/M = 0.9 




QDPa; a/M = 0.9 




TDPa; a/M = 0.9 




Figure 3. Time-averaged radial profile of Qmw (solid lines) compared to Qnt (dashed lines) and Qak (dot-dash lines) for (from left to 
right) simulations KDPg (dipole magnetic field topology), QDPa (quadrupole magnetic field topology) and TDPa (toroidal magnetic field 
field topology). The black hole spin was fixed at a/M = 0.9. Dashed lines denoted ±1 std. deviation from the mean Qmw- Simulation 
data was time-averaged over the periods given in Table [T] 



We then set S'"^ — Tj'^'m) ^^"^ calculate both this and the 
derivatives of Ufj, directly from the simulation data. This 
proced ure is similar to that used by lArmitage fc Reynolds! 
(|2003l ) in their study of the observational implications of a 
pseudo-Newtonian cylindrical disc simulation. In their pro- 
cedure they set S(^r){4,) equal to BrB^/iir as calculated from 
the simulation. 

Thus, this procedure generalizes the AK method in sev- 
eral ways. It allows for non-zero Ur (but not non-zero ug), 
and it permits a smooth extension of the estimated dissi- 
pation across the ISCO. However, it does so at the cost of 
some internal inconsistency because T^'^'m) ^'-'^ orthogonal 
to M^, a condition required by the model equations. Fortu- 
nately, this inconsistency may be of limited magnitude in 



the region of greatest inte r est, th e marginally stable region. 
As shown bv lKroUk eral] (|2005h . \\b\\' vTu^ integrated over 
the volume occupied by bound material is in most cases at 
least an order of magnitude less than —b^b^ integrated over 
the same region, and only deep inside the plunging region 
or, at very high spin {a/M = 0.998) inside r ~ 2M, does 
the advected magnetic energy become close to being com- 
petitive. Nowhere does it exceed the conventional magnetic 
stress. 

It should be noted, however, that by excluding those re- 
gions containing bound matter but outside one density scale 
height from the plane we are excluding any dissipation that 
might contribute to "coronal" emission. We do so in order to 
focus attention on the radiation coming from the disc body. 
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Figure 4. As in Figure[T]for tlie time-averaged fractional increase 
in dissipation rate, (ADmw) derived from Qmw- 



Figure 5. The fraction of the dissipation Qmw that occurs out- 
side the ISCO, r > r-rriB: as derived from time-averaged simulation 
data. The symbols have the same meaning as in Figure [T] 



the component most likely to be thermalized. Nonetheless, it 
is possible that the "coronal" portion we neglect here could 
be significant in the total energy budget of the system. 

In the simulations, only the portion of the flow lying in- 
side of r = 20M is accreting in a quasi-steady state manner; 
outside this point matter flows out as it absorbs angular mo- 
mentum transported outward from the inner disc. For this 
reason, it is inappropriate to include regions with r ^ 20M 
in the computation of Qmw- However, standard disc models 
dissipate 20-60% (for a/M = 0.998-0.0) of the total heat in 
the region r > 20M; this is clearly non- negligible. We as- 
sume, therefore, that the total energy dissipated in the sim- 
ulated disc outside of r = 20A'I matches that specifled by 
Qak- 

Examination of the radial profile of Qmw reveals that in 
the disc body it is often offset by a small amount, sometimes 
positive, sometimes negative, relative to Qak- We there- 
fore renormalise Qmw so that it exactly matches Qak at 
r = lOM in each case. The shift is generally only a few 
to ten percent, and it ensures that both models predict the 
nearly the same dissipation in the outer disc. The error it 
induces in the dissipation rate in the inner accretion flow is 
small compared to the other uncertainties of our method. 
Also note that this re-normalization is irrelevant to the lo- 
cation of the dissipation edge. For this quantity, all that is 
important is that the dissipation profile has a continuous 
extension beyond the r — 20M surface. Outside this radius, 
Qak follows an approximately scaling, similar to that 
of Qmw inside of r = 20Af and as such is an appropriate 
choice for locating the radiation edge. 

Figures [2] and [3] show the radial profile of Qmw com- 
pared to Qnt and Qak for (respectively) a variety of black 
hole spins and magnetic field topologies. Outside of the 
ISCO, Qmw, Qnt, and Qak are all very similar. However, 
as the ISCO is approached from the outside, although Qmw 
and Qak track each other well, both separate from Qnt 
because the artificial boundary condition imposed in the 
NT model forces the stress to die there, whereas the sim- 
ulations show that it in fact cont inues to rise inward (as 
also seen in the data presented in iKrolik et al] (|2005l ) and 
iBeckwith et al.l l|2008h '). Inside the ISCO, only Qmw is de- 



fined; in all cases, it rises steeply with diminishing radius 
throughout the plunging region. 

There is generally a close correspondence between Qmw 
and Qak for all radii where both are defined. This is not 
surprising, given that both models share a similar origin 
and the AK model is calibrated by the value of the stress 
at the ISCO found in the same simulations used to deter- 
mine Qmw - In addition, we further constrain Qmw to match 
Qak at r = lOM. There are some notable differences in the 
curves, however. In the high-spin KDE model, for example, 
Qak is consistently greater than Qmw by factors of sev- 
eral for almost all radii r^s < r < 4M. There are several 
possible reasons for a lack of agreement. First there is the 
aforementioned possibility that the simulations are not in a 
true steady state. Indeed, in most cases, even where the Qak 
curve differs noticeably from Qmw, it is nonetheless within 
the lines indicating the one standard deviation fluctuations 
in the simulation data. The lack of radial dependence in the 
time-averaged accretion rate argues that the simulations are, 
however, close to steady state in an averaged sense. Second, 
there is the lack of orthogonality between and T^^j^j^ 
(as described earlier in this section). However, comparison 
of Qmw and Qak calculated by replacing with b^h" 

(so that u^S^" = is satisfied) reveals similar offsets be- 
tween the different dissipation profiles. Third, the velocities 
in the simulations are not purely Keplerian. There is a sig- 
nificant net radial component that becomes larger near and 
inside the ISCO. In particular, non-zero radial velocity in- 
troduces a new term in the definition of g'' not present in 
the AK or NT models (see eqn. I19|l . Lastly, we have iden- 
tified S^" with the Maxwell stress, neglecting the Reynolds 
stress due to correlated radial and angular velocity fluctua- 
tions, but it can also contribute to the total stress. Although 
the magnitude of this contribution to the Reynolds stress is 
difficult to determine within a global simulation, local shear- 
ing box simulations have shown that it is typically near or 
below one third the Maxwell stress in the main disc body 
(|Hawlev et al.lll995l l. 
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Figure 6. The time-averaged fractional increase in dissipation 
derived from Qmw where the region under consideration is lim- 
ited to that lying outside the ISCO, r > rms. The symbols are 
as defined in Figure [T] 



Figure 7. The location of the dissipation edge in relation to the 
location of the ISCO, rms, derived from Qnt (black crosses), 
Qak (green symbols) and Qmw (red symbols). The symbols cor- 
respond to different simulations as given in Figure [Tl 



Figure |4] shows the fractional increase in dissipa- 
tion rate, (ADmw) derived from the Maxwell stress, us- 
ing a definition analogous to eqn. 1131 where Dmw = 
/°° / QtAwdx^'^^ dx^^^ dx^^\ and dx^^^ are the coordinate el- 
ements measured by an observer comoving with the fluid 
with ti° — (it*, tt*^, 0, it*^). We find values of (ADmw) from 
10 to 200%, with most in the range 50-100%. In contrast 
with (ADak), there is little dependence on black hole spin. 
For topologies initialised with non-dipolar magnetic fields, 
(ADmw) is smaller than for those with initially dipolar field. 

In both the AK and the NT disc models all of the dis- 
sipation in the disc occurs at or outside of the ISCO. In the 
simulations, however, the plunging region makes a contri- 
bution to the total dissipation, as can clearly be seen from 
Figures [2] and |3] Figure [5] shows the fraction of the total 
dissipation for each of the simulations that comes from out- 
side the ISCO. The disc above the ISCO contributes from 
60-90% of the total energy dissipation. Comparing this fig- 
ure with Figure |3] shows that, in general, the greater the 
fractional increase in total dissipation, the more significant 
the plunging region is compared to the rest of the disc. As 
a group the dipole simulations tend to have larger stress in 
the plunging region. 

In Figure[6]we return to the quantity (ADmw), but now 
include only that portion of the disc outside of the ISCO. 
Removing the contribution from the plunging region reduces 
(ADmw) to be on the order of tens of percent. There is 
perhaps a slight trend towards larger numbers with larger 
spin. Removal of the contribution from the plunging region 
reduces the contrast between the different magnetic field 
topologies. This is mainly because in the dipole simulations 
the profile Qmw steepens significantly below the ISCO (as 
can be seen from Figure [3|. 



3.3 The Dissipation Edge 

A way to characterize the relative distribution of Q in all 
three of the dissipation models is to compute the "dissipa- 
tion edge" for each model, defined to be the radius outside of 
which 95% of the total energy dissipated within the disc has 
been deposited. Figure [7] displays the location of the dissipa- 
tion edge for Qnt, Qak and Qmw in terms of the radius of 
the ISCO, Tma- The dissipation edge for the NT model lies 
between 1.2 (for a/M = 0.998) and LSr^s (for a/M = 0.0). 
Since the stress goes to zero at the ISCO in this model and 
reaches a maximum value at some location well outside of 
this point, the majority of the dissipation also occurs well 
outside the ISCO. Consequently, radiation from these discs 
does not probe (except by photon propagation) the deepest 
regions of the black hole's gravitational potential. 

A non-zero stress at the ISCO changes the picture. In 
Figure [7] we see that the dissipation edge moves closer to 
the ISCO than in the standard model, although in the AK 
model it obviously must remain at a location > Vms- The 
AK dissipation edge is located between ~ LOr^s for a/M = 
0.998 and ~ 1.5rms for a/M — 0.0. As would be expected 
from Figure [1] changes in magnetic field topology have little 
influence on the location of the dissipation edge derived from 
Qak. 

Finally, we turn to the dissipation edge derived from 
Qmw- With stress permitted inside of the ISCO, the dissi- 
pation edge consistently lies either at or inside rms ■ There is 
a visible dependence on black hole spin. The dissipation edge 
ranges between ~ 0.4 and l.Orms as a/M increases from 0.0 
to 0.998. There is little apparent dependence on magnetic 
field topology. Some of the dependence on spin may be an 
artifact of the simulations. The GRMHD code uses Boyer- 
Lindquist coordinates, and there are fewer zones inside of 
rma for high values of a/M. As we shall see in the next 
section, however, even when the dissipation edge is deep in- 
side the plunging region, the radiation edge is generally well 
outside, as few photons can escape from so near the horizon 
when the black hole spin is this rapid. 
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Figure 8. The time-averagod fractional increase in luminosity, (ALak) derived from Qak for a distant observer located at do = 5° (top 
left panel), 30° (top centre panel), 45° (top centre panel), 60° (bottom left panel), 75° (bottom centre panel) and 85° (bottom right 
panel). Symbols correspond to different simulations as defined in Figure [T] 



4 OBSERVED QUANTITIES 



4.1 Fractional Increase in Luminosity 



In the preceding sections, we have seen how stress near and 
inside the ISCO can significantly increase the total dissi- 
pation in the accretion flow and alter the location of the 
dissipation edge. In the NT model, all but 5% of the total 
energy is dissipated outside of 1.2-1.8rms. By allowing for 
non-zero stress at the ISCO, the AK model moves the ra- 
dial coordinate of the dissipation edge inward by factors of 
~ 0.8 compared to the NT model. Relaxing the assumption 
that dissipation terminates at the ISCO causes the dissipa- 
tion edge to move inside the ISCO. From an observational 
standpoint, these distinctions can be important. They affect 
the degree to which a given accretion flow can probe a black 
hole's gravitational potential, and also determine such ob- 
servables as the peak in the disc thermal spectrum and the 
width of the Fe Ka line profile. 

While the fluid frame dissipation is important, the real 
test is whether these differences carry over into the rest 
frame of a distant observer. In this section we examine this 
issue using the photon transfer function, which incorporates 
effects due to gravitational and orbital Doppler shifts, along 
with the influence of light-bending and gravitational lensing. 

We will examine the properties of these discs in two 
ways; first, when observed from different inclination angles, 
and second by averaging over solid angle to obtain a to- 
tal value. The solid angle averaged value of some observed 
quantity nio is defined as 

/ / mo(^o) sin 9oddod(j} 



{M} 



(20) 



J J sin 9od0od(fi 

The solid angle average of a given quantity is not itself ob- 
servable, and for many potential observables the dependence 
on viewing angle is considerable. Nevertheless, the angle av- 
erage provides a simple way to summarise the impact of the 
effects we are studying. 



For a given dissipation model, we can compute the total 
luminosity carried to a distant observer using the observed 
flux, eqn. [S] and integrating: 



/ / 5'[<.rf (a,/3)]Q[<urf(«./?)]x 
5{Eo - gE^)dadP; 



(21) 



Q corresponds to Qak, Qnt or Qmw as appropriate. We can 
then compute the fractional enhancement of the luminosity 
over the NT model using 



, , , iAK.MW — iNT , 

^A_Lak,mw/ = \ /. 



(22) 



The results for the AK model are plotted in Figure [8] 
for a distant observer at inclinations ranging from 5-85° . In 
i]3.1l we found (ADak) was strongly correlated with black 
hole spin, ranging from 10-500% as a/A4 approached one. 
Figure |8] shows that the corresponding fractional increases 
in luminosity are more modest, ranging between 5-40% for 
6o ^ 45° (low inclinations) to substantially above 40% for 
the half of solid angle with 6o ^ 60° (high inclinations). At 
low inclinations, (ALak) exhibits little dependence on black 
hole spin. At high inclinations, (ALak) exhibits a similar (if 
slightly weaker) dependence on black hole spin to that seen 
in (ADak). 

Figure [5] gives (ALmw) for the same range of inclina- 
tions. In i]3.2l (see Fig. |4]), we found {ADmw) was generally 
within a factor of 2 of 100%, with little consistent depen- 
dence on black hole spin. Figure [5] shows a different picture. 
At all inclinations, (ALmw) at the highest spin is a few tens 
of percent, varying by ~ 2 for the different magnetic topolo- 
gies. At all inclinations, there is also a rise in (ALmw) with 
diminishing spin, but the slope of this rise increases sharply 
with increasing inclination angle. As a result, for angles of 
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Figure 9. As in Figurc|8]for the time-avcraged fractional increase in luminosity, (ALmw) derived from Qmw- 



75° and more, (ALmw) > 100% for all spins a/M < 0.9. 
There is a weak tendency for a/M = 0.5 to yield the great- 
est fractional enhancement in luminosity, but it is not by a 
wide margin over the spinless case. 

In i)3.2l we found that the fraction of the total energy 
dissipated outside the ISCO ranged from ^ 60% to ^ 80% 
for non-dipolar magnetic field topologies. Figure[TD]plots the 
fraction of the total luminosity that is received from outside 
of Tms for Qmw- This fraction is near 100% for low incli- 
nations and a/AI > 0.9. As the spin of the hole is reduced 
below a/AI — 0.9, the fraction falls to > 90% as more pho- 
tons escape from inside the ISCO. At higher inclinations 
{6o ^ 60°), the contribution from the plunging region in- 
creases for all but the most rapidly rotating black holes, 
in which > 80% of the total observed luminosity originates 
outside of the ISCO. Thus, the contribution of the plunging 
region is greatest for low to moderate black hole spins: at 
high inclination the plunging region contributes 20-60% of 
the total observed luminosity of the disc for this range of 
spins. 

The difference in solid angle averaged luminosity com- 
pared to the standard NT value is shown in Figure [TT] Both 
{ALak} and {ALmw} are 10-100%. {ALak} generally in- 
creases with increasing a/M . There is considerable variation 
between simulations with different magnetic field topologies 
but with the same black hole spin. The AK model is sensi- 
tive to the precise value of the stress used at the ISCO, but 
even for the weakest stress values used here there are still 
significant increases in total luminosity. 

In contrast, {ALmw} decreases with increasing a/M 
and is largest for a/M = 0.5 where {ALmw} ~ 100% 
In fact, {ALmw} is larger than {ALak} at low spin, but 
smaller at high spin, consistent with results obtained in the 
fluid frame. Two effects are important here. First, the plung- 
ing region is larger for low spin holes, so emission from that 
region (absent by assumption from the AK models) is more 
significant for those cases. Second, the stress tends to in- 



crease sharply at the ISCO for high spin holes, boosting the 
emission in the AK model. As we have already remarked in 
other contexts, different magnetic topologies produce signif- 
icant variation among the simulations performed at a given 
black hole spin. 

4.2 The Radiation Edge 

The standard NT model is highly constrained. The assump- 
tion that stress goes to zero at the ISCO, the location of 
which is determined by the mass and spin of the hole, sets 
the location of the radiation edge. Removing this boundary 
condition means that the radiation edge can move, but by 
how much? As emphasised in the Introduction, one prin- 
cipal goal of this analysis is to locate the radiation edge 
for models that have nonzero stresses at the ISCO. While 
one might define the radiation edge as that point outside of 
which 100% of the light is emitted, such a definition fails 
to distinguish between the AK and NT models because for 
both this point is the ISCO. To bring out the distinctions in 
the distribution of the emission near the inner edge we de- 
fine the radiation edge as the radius outside of which 95% of 
the light that reaches distant observers is emitted. To further 
develop these contrasts, we also examine the cumulative dis- 
tribution function for the fraction of the solid angle-averaged 
luminosity reaching infinity from outside a given radius. 

In Figure [12] we plot the radiation edge for each of the 
dissipation models and for the same range of 9o shown in 
Figures 151 4101 The first thing to note is that the transport 
of photons to infinity introduces a significant separation be- 
tween the radiation edge and the previously computed dis- 
sipation edge. 

For the NT model at angles 9o ^ 45° , the radiation edge 
lies between 1.7-2.8rms, increasing with increasing a/M, 
whilst for 9o ^ 60°, the radiation edge becomes approx- 
imately independent of black hole spin and can be found 
between 1.3 and l.lrms - Comparing the location of the radi- 
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Figure 10. As in Figurc|8]for the fraction of the total luminosity emitted outside of the ISCO derived from Qmw- 



ation edge to the dissipation edge reveals that, for Oo ^ 45°, 
the radiation edge is either at (for spins a/M < 0.9) or out- 
side (for a/M ^ 0.9) the location of the dissipation edge. 
When the inclination is greater {9o ^ 60°), the radiation 
edge lies either inside [a/M < 0.9) or near {a/M > 0.9) the 
location of the dissipation edge. 

Similar behaviour is found for the AK radiation edge. 
The principal difference is that, in this case, the radiation 
edge lies inside that of the NT by a fraction that remains ap- 
proximately constant as 9o is varied, with variations only on 
the ~ 5% level for a given simulation. For example, for simu- 
lation KDOc, the AK radiation edge lies between 0.83 — 0.80 
of the NT radiation edge as 6o varies between 5° and 85°, 
whilst for simulation KDPg, this fraction ranges between 
0.80 — 0.75 over this same range of 9o- The larger the black 
hole spin, the greater the difference between the AK and 
the NT models. In other words, the relative visibility of the 
region near the ISCO remains the same function of spin and 
inclination angle for both models; the AK model simply has 
greater luminosity near the ISCO by an amount that in- 
creases with black hole spin. 

We next consider the radiation edge derived from Qmw . 
In i j3.3l we found that the location of its dissipation edge 
was inside that of the AK or NT model for all spins and 
field topologies. By contrast, the location of the radiation 
edge shows a wider range of behaviours. At any given in- 
clination, the radiation edge moves to larger r/r„is as the 
spin increases. For a given spin, the radiation edge moves 
inward as the inclination grows. Specifically, for discs that 
are viewed almost face-on {9o =5°), the radiation edge lies 
between l-Srms , increasing with increasing black hole spin. 
At this inclination, the radiation edge lies approximately a 
factor of three outside the location of the dissipation edge, 
independent of black hole spin. At moderate inclinations 
(30° ^ 6^0 ^ 45°), the radiation edge moves inwards towards 
the ISCO, ranging between 0.5 and 2.5r-ma, and again in- 
creasing with increasing black hole spin. For Oo < 45° and 



a/M > 0.9, the radiation edge lies either at or outside that 
derived from Qnt and Qak, whilst for a/M < 0.9 the ra- 
diation edge lies either at or inside that derived from Qnt 
and Qak. 

These trends in the location of the radiation edge, and 
in particular, the way it is offset from the dissipation edge 
f H3.3|l are driven by the way photons travel through the rel- 
ativistic potential. When the disc is viewed face-on, gravita- 
tional redshift dominates, reducing both the energy of pho- 
tons emitted close to the hole and the apparent rate at which 
they are released; consequently, the radiation edge tends 
to move outward as the inclination angle becomes smaller. 
As the viewing direction moves closer to edge-on, Doppler 
boosting of the approaching side of the disc increases the 
energy of photons radiated there and their rate of emission, 
bringing the radiation edge inward. In the NT and AK mod- 
els, all velocities are azimuthal, so the peak approach veloc- 
ity occurs for matter passing through the sky plane when our 
view is nearly in the disc equatorial plane; however, in the 
MW model, the inward radial speed can be great enough, 
especially in the plunging region, for Doppler boosting also 
to enhance the luminosity of matter on the far side of the 
black hole. In all models, emission from the far side of the 
hole is also strengthened by gravitational lensing when the 
line of sight is near the disc plane. Thus, the radiation edge 
moves inward as the observer approaches the disc plane. 

Trends in the position of the radiation edge as a function 
of black hole spin at fixed inclination angle are the result of 
a different trade-off. As the spin increases, the ISCO moves 
inward and all relativistic effects are strengthened. Those 
described in the previous paragraph tend, on balance, to 
make the radiation brighter as the inclination angle grows. 
However, as the spin increases, the fraction of all photons 
captured by the black hole also increases, likewise because 
the ISCO moves inward. In terms of how far inside the ISCO 
the radiation edge falls, the latter effect is the strongest: the 
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Figure 11. Solid angle averaged {ALak} (left panel) and 

ratio of the radial coordinate of the radiation edge to rma is 
least for low-spin black holes viewed nearly edge-on. 

One surprising result demands special discussion: at the 
highest spin (a/M = 0.998) and small inclination angle 
{60 ^ 30°), the MW radiation edge either coincides with 
or lies a little bit outside the NT edge, and both the MW 
and NT edges are well outside the AK edge. At such high 
spin, photons escaping to infinity in the polar direction are 
severely redshifted if their point of origin is near or inside 
the ISCO. For this reason, all three radiation edges in this 
regime are at radii several times r^a ■ As we have already dis- 
cussed, however, when the spin is this rapid, the AK model 
predicts a rather larger dissipation rate in this range of radii 
than is predicted by either MW or NT. The disparity be- 
tween the AK radiation edge on the one hand and the MW 
and NT edges on the other follows directly from this con- 
trast. In addition, the radiation edge for one high-spin sim- 
ulation (KDEa) appears to fall at particularly large radius, 
although this may well be an artefact o f the lack of inflow 
equilibrium in this particular simulation (jKrolik et al.ll2005l ) 

The radiation edge associated with the solid angle- 
averaged luminosity is shown in Figure 1131 For Qnt, it 
lies between 1.5-1.8rms and slowly increases with increasing 
black hole spin. For all models as the spin increases, photons 
radiated deeper in the potential are subject to larger grav- 
itational redshifts and are more likely to be captured. For 
Qak, the radiation edge moves inwards to 1.3-1.6rms, with 
the smallest values occurring for a/M — 0.9 and increasing 
as one moves away from this spin in either direction. For this 
model the enhanced dissipation at the ISCO compensates 
somewhat for the greater likelihood for photon capture with 
increased black hole spin. For Qmw, the radiation edge lies 
between 0.5-1.7r„is, increasing with increasing black hole 
spin. The enhanced dissipation within the plunging region 
boosts the luminosity (and reduces the radius of the radia- 
tion edge) for low spin models, but has little effect for high 
spins. 

Regardless of spin the solid angle- averaged radiation 
edge derived from the simulation data lies within that of the 
NT model; enhanced stress always moves this point inward. 
Comparing Qak and Qmw we find a more complex picture. 
The radiation edge in the AK model must, by assumption. 
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{ALmw} (right panel). Symbols are as described in Fig. [T] 

lie outside the ISCO, but for the models with stress in the 
plunging region this need not be the case. Indeed, for low 
spin holes the radiation edge can be inside the ISCO, but 
the dissipation in the plunging region becomes less impor- 
tant as black hole spin increases. What matters most is the 
dissipation level outside of the ISCO. For the cases consid- 
ered here, the stress rises rapidly at the ISCO for high-spin 
models and the semi-analytic AK formula predicts greater 
dissipation outside the ISCO than is implied by the stress 
levels seen in the simulations. Hence, for a/M > 0.9 the 
solid angle-averaged radiation edge derived from Qmw lies 
outside that derived from Qak- This clearly illustrates that 
the location of the radiation edge can be very sensitive to 
the dissipation levels near the ISCO. 

We have defined the radiation edge as the the radius 
outside of which 95% of the total observed luminosity is 
emitted. This choice was made to provide a simple mea- 
sure by which to contrast the different luminosity profiles 
associated with each of the dissipation models. In Figure 
1141 we plot the solid-angle averaged fractional cumulative 
distribution of observed luminosity, {F(r/r„is)} for each 
of the dissipation models and individual datasets. We see 
that for 0.2 ^ -F ^ 0.95, all these curves may be fit- 
ted to a reasonable degree of approximation by the form 
{F{r/r,ns}) oc exp[—{r/rms)/r,]. Thus, there is a charac- 
teristic radial scale for the luminosity profile that is always 
related to, but not necessarily identical to, r^s, and it may 
be adequately parameterized by choosing a fiducial level of 
F. 

The plots also provide a quantitative sense of the poten- 
tial importance of nonzero stress at the ISCO. First, note 
that the luminosity profile in the NT model is practically 
independent of spin (as can also be inferred from the solid- 
angle averaged radiation edge shown in Figure I13p . This 
somewhat surprising result can be understood in terms of 
the dependence of the radiation edge on 60 ■ At small incli- 
nation angle (face-on views) , the radiation edge derived from 
the NT model for slowly spinning black holes lies inside that 
of the rapidly spinning holes (relative to the location of the 
ISCO), whilst at large inclination angles the reverse is true. 
When this data is averaged over solid angle, these changes 
cancel and so the position of the radiation edge (and hence 
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Figure 12. Timo-avcraged "radiation edge" seen by distant observers in in six different directions (arranged as in Fig.|8] do = 5° (top left 
panel), 30° (top centre panel), 45° (top centre panel), 60° (bottom left panel), 75° (bottom centre panel) and 85° (bottom right panel). 
The different models are coded by differently coloured symbols: Qnt (black), Qak (green) and Qyiw (red) . Different symbol shapes 
denote differ ent simulations, as for Fi gure [T] The dipolar field simulations described bv lDe Villiers et al.l 1120031 ) are shown by stars, those 
described bv lHawlev fc KrolikI l|2006h by diamonds and the new high spin cases computed for this paper by triangles. Simulations with 
quadrupolar topologies are shown by squares. The toroidal field simulation (a/A/ = 0.9) is denoted by a cross. 



the fractional luminosity distribution) relative to r,„s is in- 
dependent of a/M. The NT model is so tightly constrained 
by its assumptions that when the dependencies due to ob- 
serving angle are removed, there remains almost no contrast 
other the location of the ISCO which is set by the spin of 
the hole. 



4.3 Characteristic Temperature 

We have seen how the apparent size of the disc varies 
with spin, inclination and dissipation model. However, ob- 
servations do not directly measure the radiation edge. 
Rather, the concept of a radiation edge is incorporated 
into the interpretation of the L = A(a/M,6)T^ relation- 
ship (|Gierliiiski fc Don^ 12004 ). We have already seen how 
changes in the dissipation function due to non-zero stress 
at and inside the ISCO have the potential to change the 
apparent area and luminosity of the disc. In this section, 
we attempt to gauge the impact of these changes on the 
characteristic temperature of the accretion flow, Tchar, as it 
might be inferred from continuum fitting to the soft com- 
ponent of a spectrum. Our analysis must necessarily be 
simple. We define Tchar as the maximum (observed) black- 
body temperature found anywhere in the disc. It is de- 
termined by defining an effective temperature in the fiuid 
frame: Te//(r) = [Q{r) / a]^^'^ and then transforming this 
temperature to t he rest frame of a distant observer via 
To{r) = gTeffir) (|Cunninghamlll975l '). Then, since the char- 
acteristic temperature of the associated blackbody spec- 
trum will be very close (but not identical to) the maxi- 
mum black body temperature, we make the approximation 
Tchar ~ max(ro). We perform the same procedure for all 
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Figure 13. Solid angle averaged {r^dge/fn 
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three dissipation models, ignoring color temperature correc- 
tions and all observational effects except those due to incli- 
nation angle. 

Although the geometric size of the disc in the AK model 
is the same as in the NT model (the curves are still con- 
strained to reach 100% at the ISCO), the middle panel of 
Figure [T4l shows how the addition of a finite amount of stress 
at the ISCO breaks the NT degeneracy. The most signif- 
icant differences come, naturally, from including contribu- 
tions from inside the ISCO. The curves for the Qmw models 
(right hand panel) are well separated for different spins. For 
these models, the relative importance of the plunging region 
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Figure 14. Solid-angle averaged fractional distribution of observed luminosity, {Fir /rms)} derived from Qnt (left-hand panel), Qak 
(centre panel) and Qmw (right-hand panel). 



is a function of both the amphtude of the stress there and 
the spin of the hole. The farther the ISCO from the hori- 
zon, the greater the chance for radiation within the plunging 
region to escape. 

Tchar is plotted in Figure [15] for the different dissipa- 
tion functions and observer inclinations shown in Figure [T^ 
Each of the accretion flows has been scaled such that its ac- 
cretion rate would produce a flow with luminosity O.lLedd. 
Several features are immediately apparent. At all inclina- 
tions Tj^Jf^, the characteristic temperature associated with 
Qnt increases with increasing a/M. This well known re- 
sult is the foundation of attempts to measure black hole 
spin with spectral fitting techniques. T^^^ shows a simi- 
lar trend, with an increase over Tj/J^r by a factor 1.1 — 1.4, 
an amo unt that is comparab le to typical color corrections 
(see e.g. lDone fc Davi3l2008h . The AK model has the same 
(fixed) inner boundary, but has enhanced luminosity for the 
same accretion rate. The effect is most significant for large 
inclination angles. 

Tch^ exhibits a different behaviour. At low spins, T^h^ 
is significantly greater than rj^an whilst at high spins the 
contrast is reduced. Overall, this has the effect of making 
T^jl^ (a/M) almost constant for a given inclination. Over- 
all, T^Z /TSZt varies between a factor of 1.8 - 2.3 for 
accretion flows accreting at the same luminosity, with the 
greatest increases occurring for slowly spinning black holes 
viewed nearly edge on. Because energy extraction and radi- 
ation occur within the plunging region, the distinctions be- 
tween holes with different spins is greatly reduced. In other 
words, the effective inner boundary for all discs lies close to 
the horizon. 

These effects are summarised in Figure 1161 where we 
plot the characteristic temperature associated with the 
solid-angle averaged radiation edge for a black hole accret- 
ing at O.ll/edd- Both the NT and the AK model share the 
ISCO as their inner boundary, and as a result both Tj^or 
and T^^^ increase by around a factor of two over the range 



oi a/M shown here (both increasing with increasing a/M). 
On the other hand T^hZ is approximately constant over the 
same range of spin and is greater than by a factor that 

falls from ~ 2.3 to ~ 1.2 as a/M increases from to 0.998. 
Unlike the radiation edge, the contrast between Tjf^ and 
both and T^'^^ always has the same sense: the MW 

model always predicts a higher characteristic temperature. 

4.4 Impact on Measurements of a/M 

In the previous section, we examined the consequences of 
stress at the ISCO on the characteristic temperature of the 
accretion flow and the dependence of this quantity on both 
black hole spin and inclination. We now examine this ques- 
tion from the opposite direction, i.e., supposing one has a 
measurement of the characteristic temperature of an accre- 
tion disc in a particular system (along with the black hole 
mass, inclination of the binary orbit and disc luminosity), 
how much uncertainty in the determination of a/M is in- 
duced by uncertainty about the correct dissipation profile 
in the inner disc? 

To address this question, we use the data shown in Fig- 
ure [15] to construct Figure [17] in which four polygons, one 
for each of four different inclinations, illustrates the range 
of Tchar and a/M consistent with the different disc mod- 
els. To produce this figure, we first fixed the inclination, 
spin, and black hole mass, and then collected the values of 
Tchar predicted on the basis of any of our three candidate 
models (NT, AK, and MW) from any of the several simu- 
lations we conducted with that particular spin parameter. 
The left-hand edge of each region traces the prediction of the 
NT model because it always gives the lowest temperature. 
The horizontal width of each region is defined by the range 
of temperatures predicted by the complete complement of 
models and simulations. Because we have conducted simu- 
lations at seven different spins, the right-hand edge is defined 
at seven points. 
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Figure 15. Timc-averagcd characteristic black body temperature for a IOM0 black liolc accreting at O.lLe^d measured by distant 
observers in in six different directions (arranged as in Fig. 1121 with the same key as that figure). 



Several qualitative features emerge from study of this 
figure. First, for each inclination, the maximum tempera- 
ture change predicted by the NT model as the spin runs 
from a/M = to a/M = 0.998 (generally about a factor 
of 2) is always less than the maximum temperature change 
at zero spin due to model-dependent uncertainty in the dis- 
sipation profile (typically a factor of 3). Thus, if we knew 
the inclination angle and mass for a given black hole, but 
not its spin paramater, our uncertainty in predicting Tchar 
has a greater contribution from our uncertainty about disc 
physics than it does from our ignorance of its value of a/M . 

Second, suppose that we know the inclination angle of 
a particular black hole and obtain Tchar from X-ray spectra. 
The uncertainty in a/M due to uncertainty in the dissipa- 
tion profile may be estimated by imagining a vertical line 
at that value of Tchar running through the region for the 
appropriate inclination angle. The uncertainty is then given 
by the difference in spins between the points where the line 
intersects the boundary of the region. The size of this un- 
certainty tends to be larger when our view is more nearly 
edge-on. However, for almost any value of Tchar that is con- 
sistent with any of the area of a region, the magnitude of 
the uncertainty in 1 — a/M is generally at least an order of 
magnitude. 

Still another way to look at this diagram is to ask, "If 
we estimate spin by using the NT model, by how much may 
it be wrong if a different model is truer to the disc's ra- 
diation profile?" The answer is the span of a vertical line 
whose lowest point is the NT prediction and runs to the 
upper boundary of the associated region. Because the NT 
model always gives the lowest prediction for Tchar, the spin 
inferred from it is always the greatest possible spin, and the 
error bar always stretches toward lower rotation rates. Given 
the characteristic curvature of the permitted regions in this 
plane, the magnitude of the possible error in 1 — a/M tends 
to be larger (in logarithmic terms) when the NT inference 



of a/AI is closer to unity, but can still be substantial even 
when a/M < 0.9. Note, however, that we only present data 
for a/M > 0, so we cannot provide an explicit lower bound 
on the range of a/M , beyond stating that a/M = in the 
magnetized case is always consistent with the full range of 
a/M in the NT case for fixed 60- 



5 SUMMARY, DISCUSSION AND 
CONCLUSIONS 

That magnetic forces might cause substantial stress at the 
ISCO was fores een very shortly after the invention of the 
standard model (|Page fc Thorndl 19741 ). This possibility now 
appears to be an immediate corollary of the well-established 
result that MHD stresses account for most of the angular 
momentum transport in the bodies of accretion discs. In- 
deed, such stresses were seen in the first generation of three 
dimensional GRMHD disc simulations. The goal of this pa- 
per has been to begin the linkage of these numerical MHD 
simulations to the observable properties of accreting black 
holes even before the simulations are fully equipped to make 
predictions about how these systems radiate. To do so, we 
have followed a path of cautious extrapolation from older 
methods. We first used simulation data to fix the single 
parameter of a model (called AK here) that changes the 
previous standard (the Novikov-Thorne model) only by ad- 
mitting the possibility of a non-zero stress at the ISCO. 
Because the AK model is defined in a way that prevents 
its extrapolation within the ISCO, we used the formal- 
ism underlying both it and the NT model (i.e., vertically- 
integrating and azimuthally- and time-averaging the equa- 
tion of momentum-energy conservation under the assump- 
tion that the four-velocity and the stress tensor are orthogo- 
nal) to create an expression for the dissipation (called Qmw 
here) valid both inside and outside the ISCO. 
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Figure 16. Plot showing solid angle averaged {T^har} • Symbols 
are as in Figure [7] 



Happily, in the body of the accretion disc well out- 
side the ISCO, both the AK and the MW method agree 
fairly well with NT more or less independent of black hole 
spin and magnetic field topology, although the irregularity 
in the curves of Figure [5] reminds one that the simulations 
are dynamic and time-varying, and that the 26 samples of 
simulation data we used define only somewhat imperfectly 
the long-term time-average. In addition, with the excep- 
tion of the extreme high-spin example, where the AK and 
MW methods depart from NT just outside the ISCO, they 
do so together. This is noteworthy because, although they 
are based on closely-related formalisms, they are not iden- 
tical: perhaps their most significant contrast is that the AK 
method assumes that Ur = 0, while the MW method does 
not. Lastly, inside the ISCO, where only the MW method 
is defined, it follows a smooth extrapolation from larger ra- 
dius. When the black hole spins slowly [a/M < 0.9), Qmw 
extends with hardly any change in logarithmic derivative 
with respect to radius. For higher spin, the extension grad- 
ually steepens toward smaller radius, but the next step in 
our formalism shows that this makes little difference to ob- 
served radiation; relativistic ray-tracing shows that the vol- 
ume deep inside the plunging region, particularly at high 
spin, contributes little energy to the luminosity reaching ob- 
servers at infinity. Thus, we are relatively confident that, 
despite the uncertainties involved, our estimate of the loca- 
tion of the radiation edge is comparatively insensitive to the 
exact relation between the fiow's detailed dynamical prop- 
erties and the dissipation rate. 

The dependence of the radiation edge on spin may be 
summarised succinctly: At the highest spin, there is rela- 
tively little difference between the different methods of esti- 
mating its position because the ISCO is so close to the hori- 
zon that the great majority of photons released in the plung- 
ing region never reach infinity, or if they do, are severely red- 
shifted. It moves from 2-3 times the radius of the ISCO when 
the disc is viewed face-on to almost exactly at the ISCO 
when the disc is viewed nearly edge-on. This inward move- 
ment of the radiation edge with increasing inclination an- 
gle is quite model-independent, as it stems from relativistic 
photon propagation effects: when photons from the plung- 
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Figure 17. Plot showing the region of the (Tcfm^, a/M)-plane 
consistent with a given value of 6o for a IOMq black hole accreting 
at O.lLetjd- Different colors correspond to different 6o (see key in 
bottom left corner of plot), dot-dash lines show the T^har curve 
predicted by Qjvt at a given 9o, and the region enclosed by the 
dot-dash and dashed lines of a given color show the region of the 
{Tchary a/M)-plane consistent with that value of 9o. 

ing region do reach infinity with substantial energy, it is be- 
cause they are emitted in the direction of the orbital motion. 
Probability of escape is then enhanced by a combination of 
special relativistic beaming and gravitational lensing; energy 
at infinity is enhanced by special relativistic Doppler boost- 
ing. As the black hole spin decreases, the diminishing depth 
of the potential immediately inside the ISCO makes it pro- 
gressively easier for photons to escape from that region and 
reduces the gravitational redshift they suffer when they do. 
The result is that the radiation edge moves farther inside the 
ISCO as either the spin diminishes (at fixed viewing angle) 
or the inclination angle moves toward the equatorial plane 
(at fixed spin). At its most extreme, the case of a/M = 
and 6o = 85°, rre can be ~ 0.5rms- 

Figure [10] gives additional cause to believe that these 
results are comparatively insensitive to dissipation model. 
Although the radiation edge can move well inside the ISCO 
at low spin and high inclination, most of the light received 
by distant observers is generally emitted in the region near 
and outside the ISCO. Only at the highest inclinations {9o > 
75°) and lowest spins {a/M <^ 0.5) does the contribution of 
the plunging region to the luminosity approach 50%. Thus, 
most of the light seen at infinity likely comes from a region 
where the predictions of the AK and the MW models differ 
little. 

Nonetheless, because it is also true that most of the 
light is emitted within a radius at most a few times the 
ISCO (except for the highest spin viewed more or less face- 
on), the contrast in total luminosity between the AK and 
MW models on the one hand, and the NT on the other, are 
order unity for all spins ^ 0.9. For higher spins, the effect 
may be smaller, but the uncertainties are also greater. 

These conclusions have immediate implications for a 
number of p h enom enological issues. Firstly, as suggested by 
iFalcke et al] (|2000D . it may be possible to image the near- 
est supermassive black hole, the one in Sgr A*. Because its 
accretion flow, unlike those of intrinsically brighter systems. 
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could well be radiatively inefficient, a simulation scheme that 
conserves total enersy is mo re appropriate to analysing its 
emission. iNoble et alj (|2007l ). using such a code (albeit an 
axisymmetric version) , have produced predicted images that 
illustrate several of the effects emphasised here, although in 
their work so far they have not reported quantitative de- 
scriptions of characteristic emission radii. 

Secondly, the enhanced total radiative efficiency due to 
dissipation in the marginally stable region may affect esti- 
m ates of population-mean spin parameter s (e.g., as for AGN 
bv lElvis et al]|2002l : IYu fc Tremainell2002h . Because the effi- 
ciency rises with increasing prograde spin in the NT model, 
the spin inferred by this method may overestimate the actual 
spin of accreting black holes if this enhancement is ignored. 

The additional luminosity from enhanced stress in the 
innermost part of the accretion flow could significantly alter 
the emergent spectrum. Employing a simple thermal model, 
we have found that the characteristic temperature of the 
flow increases by a factor of 1.2 -1.4 over that predicted by 
the NT model. As a consequence, the thermal peak of the 
disk spectrum (at ~ 1 keV in Galactic black holes, ~ 10 eV 
in AGN) may be pushed to somewhat higher energies. 

Several caveats must be mentioned, however, in regard 
to this prediction. First, this number supposes an emer- 
gent spectrum that is Planckian, but most estimates of the 
disc atmosphere's structure suggest that it is scattering- 
dominated, so that the color temperature of the spectrum 
is shifted upward from the effective temperature. The mag- 
nitude of this shift depends on details of the d isc's vertical 
struc ture that are n ot as yet well known (see e.g. lDavis et al.l 
I2OO5I ). Furthermore. iBlaes etliD (|2006l l (using the vertically 
stratified shearing box simulations of iHirose et al.l I2OO6I ) 
show that magnetic pressure support changes the vertical 
structure of the disk resulting in a noticeable hardening 
of the emergent disk spectrum compared to the standard 
Novikov-Thorne picture due to non-LTE effects. Second, it 
is possible that some of the enhanced dissipation will oc- 
cur where the density and optical depth are too low to 
accomplish thermalisation. Strengthening of the "coronal", 
i.e., hard X-ray, emission, rather than hardening the ther- 
mal disc spectrum would then be the likely consequence. 
Third, our treatment ignores those photons emitted deep in 
the potential that neither escape directly to infinity nor are 
captured by the black hole, bu t instead strike the disc. As 
shown bv lAgol fc KrolikI (|2000l ). this "returning radiation" 
can be a substantial fraction of all photons emitted when 
r < 5rg. Depending on their spectrum and the structure of 
the disc atmosphere where they strike, these photons may be 
either reflected (with Doppler shifts) or absorbed and their 
energy reradiated at a different (in general, lower) temper- 
ature. Quantitatively evaluating all three of these effects is 
well beyond the scope of this paper, but can be done in 
future work. 

There are also implications for attempts to determine 
black hole spin from spectral fitting. In all three models, the 
characteristic radius of emission is always near the ISCO, 
but does not, in general, coincide with it. Generally speak- 
ing, this characteristic radius is largest for the NT model, 
smaller (but still outside the ISCO) for the AK model, and 
smaller still, possibly moving into the plunging region inside 
the ISCO, for the MW model. Because the ISCO moves to 
smaller radial coordinate as a/M increases, these character- 



istic radii always become smaller for faster spin. However, 
the fractional amount by which the characteristic emission 
radius moves inward in the MW model is greatest for the 
lowest spins, so that in the end, the MW model predicts a 
relatively slow variation of radiation edge with black hole 
spin. The AK model, like the NT model, does not radiate 
from inside the ISCO, but the additional stress at and just 
outside the ISCO in this model (relative to the NT pre- 
diction) produces a systematic increase in the characteris- 
tic temperature. The magnitude of this shift in characteris- 
tic temperature rises, of course, with increasing additional 
stress. When there is emission from the plunging region, as 
in the MW model, the characteristic temperature can rise 
still higher, but the highly relativistic motions there make 
observed properties more strongly dependent on inclination 
angle. In addition, a larger fraction of the emitted photons 
can be captured by the black hole. 

When all these considerations are combined, we find 
that, for fixed black hole mass, luminosity, and inclination 
angle, the uncertainty in the characteristic temperature of 
the radiation reaching distant observers due to uncertainty 
in the dissipation profile is greater than the that due to a 
complete lack of knowledge of the black hole's spin. Clearly, 
our incomplete understanding of accretion disc physics (here 
specifically the magnitude of the stress at and inside the 
ISCO) makes it difficult to determine a black hole spin based 
on continuum model-fitting. The best one can say is that 
estimates based on the traditional Novikov-Thorne model 
can be expected to yield the most rapid spin possible, but 
the actual spin may be significantly slower. 

Our results demonstrate the potential importance of 
nonzero stresses at and inside the ISCO. But how repre- 
sentative are the specific values obtained in these simu- 
lated discs? There are two considerations: those arising from 
purely numerical effects, and those limitations arising from 
the assumptions and parameters of the model used. 

First, the results of numerical simulations can be influ- 
enced by finite resolution and the limitations of the numer- 
ical technique. All of the simulations presented in this work 
were performed at a resolution of 192 x 192 x 64 (r, 6, (f>) grid 
zones using ideal MHD and an internal energy equation. The 
equation of state and the numerical energy dissipation are 
unlikely to have a direct effect on our conclusions as Qmw 
is derived directly from the physical Maxwell stresses within 
the disc, rather than by measuring some numerical dissipa- 
tion rate. Low resolution usually causes the Maxwell stress 
to be undervalued; if so, the implications of this paper would 
be strengthed by improved resolution. Until available com- 
puter power makes better-resolved three-dimensional sim- 
ulations possible, the best way we have to test the effects 
of finite resolution is to compute axisymmetric simulations 
with higher resolution. A va r iety o f such simulations were 
presented in iBeckwith et al.l (|2008l ) with resolutions up to 
1024^ . We observe that greater resolution reduces the rate of 
numerical reconnection and improves the ability of the sim- 
ulation to maintain certain field configurations and small- 
scale field structures. Overall the amplitude of the turbulent 
Maxwell stresses in the disc remained largely unchanged as 
resolution was increased. We have also calculated W,'l! and 
Qmw, and find no significant qualitative differences from 
the results presented in this work. 
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Beyond the purely numerical issues, the value of the 
Maxwell stress at the ISCO may depend on a number of 
disc properties. In the ensemble of simulations presented 
here, the stress levels are determined in part by the ini- 
tial field topology (dipole versus quadrupolar, poloidal ver- 
sus toroidal, and the presence or absence of a net vertical 
field). Indeed, local shearing-box simulations suggest that 
the saturated field strength can increase s ubstantially when 
large -scale vertical field threads the disc (jBalbus Sz HawlevI 
1 19981 ). 

It is also possible that the saturation stress depends 
on disc thickness. To quantify disc thickness, we define the 
scale-height H as the proper height above the plane at which 
the time- and azimuthally-averaged density falls by 1/e from 
its similarly averaged value on the equatorial plane (see |3]) : 
H = J^'^^ \/geeir = 3M,9)de. We similarly define R as the 
proper (as opposed to coordinate) radial distance from the 
horizon to 3M plus the coordinate distance from the origin 
to the horizon: R — ri„ + J^^^^' \J grr{r, 6 = 7r/2)dr, with 
Tin as given in Table [1} ■ We then define the disc thickness 
as the ratio H/R. Even though r = 3M lies well within the 
plunging region for the lower spin cases, we find that there is 
a slow enough radial variation in this quantity to make H/R 
a reasonably well-defined parameter. Measured in this way, 
our discs are modestly thick, with a characteristic aspect 
ratio H/R = 0.06-0.2 at r = 3M. Most of the range in 
H/R results from the fact that the maximum pressure in 
the initial condition for these simulations varies somewhat 
between different a/M. 

It is difficult to say, however, just what sort of depen- 
dence there may be on disc thickness. There are some argu- 
ments suggesting that magnetic effects may increase with 
increasing H/R. For example, local shearing-box simula- 
tions find that the Maxwell stress is proportional to mag - 
netic pressure (jBalbus fc HawlevI [l9^ : ISano et al.l |2004 ). 
but there have not yet been any systematic studies of 
what regulates the magnetic pressure i n glob al, vertically- 
stratified discs. lAfshordi &: Paczvhskil (|2003l l have argued 
that inner disc stresses and dissipation may depend on disc 
thickne ss as well as on acc retion rate, an argument reiter- 
ated bv lShafee et al.l (|2007l ). but their arguments are framed 
in an essentially hydrodynamic context, and therefore elimi- 
nate any possibility of predicting magnetic stresses. They are 
also non-relativistic, and therefore eliminate any effects due 
to frame-dragging. There are also arguments that any de- 
pendence on H/R may be weak. In the simplest analytic or 
quasi-analytic MHD models, magnetic torques at the ISCO 
remain significant even in th e limit of a zero pressure disc 



Davis et al. II2OO6I : IShafee et al1l2006l : iMiddleton et al.ll2006l : 
Liu et al. 20081 ) all have spins in the range a/M ~ 0.1 — 0.8 



emam signmcant even m tn e nmit 01 a zero pressure gisc 
Kroliklll999l : lGammie|[l999l ). As discussed in iKrolik et al.l 



2005h . processes analogous to the Blandford-Znajek mech- 



anism can readily transport energy and angular momentum 
from rotating black holes to the accretion flow, and there 
is no particular reason to think that these processes should 
be tightly connected to the pressure in the disc. In the end, 
only direct simulations with the resolution to describe thin 
discs adequately will answer this question in a satisfactory 
way. The results of the present investigation provide yet one 
more reason why it will be important to do so. 

The final question that we address in this work is how 
the results presented here relate to current state of the art 
measurements of black hole spin via spectral fitting of the 
disk continuum. The 6 systems with the best data (see e.g. 



based on disk models that assume the stress-free inner 
boundary condition. From the perspective of the disk stress 
models these spins are upper limits. This might indicate 
that the hole is counter-rotating, but also opens the possi- 
bility that spin determinations might themselves constrain 
the physics near the ISCO. Firstly, the stress levels at the 
ISCO could be near the value assumed by the stress-free in- 
ner boundary condition; secondly, the classical relationship 
between stress and dissipation might not hold for enhanced 
magnetic stresses ear the ISCO; thirdly, the density levels at 
and inside the ISCO could be insufficient to thermalize the 
dissipated heat; fourthly, the time-scales for thermalization 
and radiation of the dissipated heat could be longer than 
the inffow time-scale. Another uncertainty which we have 
not examined is that the plane of the disk and the equato- 
rial plane of the hole could be misaligned and so the disk is 
subje ct to the Bardeen- Peterson effect (see e.g. [Fragile et al.l 
l2007t ). Understanding the roles played by the these effects 
will be crucial in providing robust estimates of black hole 
spins via spectral fitting of the disk continuum. 
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APPENDIX A: BASIS VECTORS OF THE 
LOCAL REST FRAME 

The calculation of the photon transfer functions requires 
the introduction of a set of basis vectors describing the lo- 
cal rest frame of the fluid (the "fluid fra me"). Such a tetrad 
set was presented bv lKrolik et al. who used a Gram- 

Schmidt orthonormalization procedure to construct it. Un- 
fortunately, some of the expressions given in that work were 
incorrectly transcribed. The correct version is as follows: 



-(0) 



k4.\J\- g'U + 94><t>9tt\ 

\g4"t>'^''' + gt4,u\o, 0, - {gtcj.u^ + fftttt*)] 



1 



kr h(h 



, I- t </> n ; r <t> 

/grrU U , , U, yffrr^J U 



/ 9rr 



where: 



kr = \/ |ff0</. (U'l')^ + Qrr {W)^ + {2gt4,U't' + QttU*) \ 

kg = \J {u't'f + Qrr {Wf + gee {u")^ + u* (2gt0if* + gttu') \ 
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